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Abstract. We propose a setup for which a power-law decay is predicted to be 
observable for generic and realistic conditions. The system we study is very simple: A 
quantum wave packet initially prepared in a potential well with (i) tails asymptotically 
decaying like ~ x~ 2 and (ii) an eigenvalues spectrum that shows a continuous part 
attached to the ground or equilibrium state. We analytically derive the asymptotic 
decay law from the spectral properties for generic, confined initial states. Our 
findings are supported by realistic numerical simulations for state-of-the-art expansion 
experiments with cold atoms. 
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1. Introduction 

The temporal evolution of an initially localized, quantum mechanical wave packet is of 
fundamental importance for our understanding of the distinctive properties of classical 
and quantum transport, and goes back to Schrodinger's early attempts to construct 
wave groups which behave alike classical particles [Tj. With increasing complexity of the 
potential landscape, which can be due to its topography as well its topology, a panoply of 
surprising transport phenomena emerge, from ballistic over sub-, super- or just diffusive 
[2]|3] to various types of localized transport jlH?] . These phenomena have their respective 
effective descriptions, e.g. in semiclassical |8|, mesoscopic [9] or statistical [TO] terms, and 
manifest in very distinct experimental settings, from light-matter interaction [TT| over 
fundamental (quantum) optics [12j[T3] to ultracold matter [H1H5], quantum walks [16] 
and biochemistry [T7j. On the fundamental level, however, all transport properties are 
hardwired in the spectral properties of the underlying quantum system, and control of 
the latter implies control of the former. Given the stunning experimental control over 
potential landscapes in variable dimensions, as achieved e.g. in cold matter science over 
the last decade [IMT8H2T] , we can contemplate to explore the above diversity of quantum 
transport phenomena for optimal control, by tuning the decisive spectral properties. 

In this paper we study the particularly simple, though paradigmatic case of a 
tunneling escape from a one dimensional potential well. The proper choice of the 
confining potential allows one to induce algebraic rather than exponential decay, for 
generic, confined initial states. As described in references [22}[23], the crucial ingredient 
will be an asymptotically scale invariant potential decaying to zero as ~ x~ 2 , which 
occurs naturally in systems with dipolar interactions [24] or anomalous molecular 
binding potentials [25J, and is often associated with quite counterintuitive effects (see, 
e.g., [22 | 1261429] ) . We thus provide in the next section [2] a particularly transparent 
and simple example of a scale-free relaxation process, where the equilibrium state is 
reached only at asymptotically long times j3l[30]. The here considered scenario is 
shown to be robust against unavoidable experimental modifications of the idealized 
theoretical scenario we depart from. In section [3] will argue that algebraic, rather than 
exponential decay, should be observable in modern experiments with expanding Bose- 
Einstein condensates [151 EO] in engineered potentials [IHltlHJET]. The final section 0] 
concludes the paper. 

2. Power-law relaxation 

We start our analysis by considering the two parameter family of ID, single particle 
potentials [27] 




x\ < L 
x\ > L , 



(1) 




(2) 
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and use dimensionless units in the sequel, measuring actions in units of H, and giving 
the particle unit mass. The exact analytical solution of the Schrodinger equation for the 
potential Vg is available [271131]. We will however see that, given its asymptotic scaling, 
the precise form of the confining well is not crucial for the power-law decay of the survival 
probability as our main result. The potential of equation (JTJ) has the same asymptotic 
behaviour as in references [22l[23], however, differently from the one considered in [23] 
it extends over the entire real axis, implying the existence of an equilibrium or ground 
state. 



2.1. Theoretical predictions 

Equation @ ensures that the ground state eigenfunction ip associated with the 
eigenvalue A = is continuous with a continuous first derivative. For |x| > L one 
finds tpo = A /x% , with a = a/1 + 8Vi — 1 [221127] . what ensures the square integrability 
for a > 1. Besides the discrete eigenvalue A = 0, the spectrum of Vg has a continuous 
component in the range E > 0, which is attached to the null eigenvalue. The associated 
eigenfunctions ipE can be written in terms of Bessel functions, and are normalized 
according to J dx iPe(x)iPe'(x) = S(E — E'). 

Let us now consider the particle prepared in an initial state ^q(x) = 0), given 
by a linear combination of eigenfunctions {ip (x), i/je{x)}, with real coefficients {a , a E }. 
Its time evolution will be given by 

POO 

tf(x,t) = a M%) + / dEa E ip E (x)e- iEt , (3) 
Jo 

and we are interested in the survival probability P(t) that the particle remains confined 
within the well, i.e., 

r+L 

P(t) = / dx\V(x,t)\ 2 . (4) 



J-L 

By substituting equation (j3J) we obtain: 

1 2 



P{t) = \a \ 2 h + 2a C 2 (t)+C 3 (t), (5) 



where 



CO 



C 2 (t)= I dEa E h{E) cos{Et) (6) 



roo roc 

/ dE dE'a E a E 'h(E,E') cos[(£ - E')t] 
Jo Jo 



/+L r+L 
dx\ip (x)\ 2 , I 2 (E) = dx^ (x)^ E (x), 
L J-L 

/+L 
dx4j E {x)^E>{x). 

For the ID potential of equation (CQ), the integrals 1%, I 2 and /3 can be solved analytically, 
while this is in general not the case for the time-dependent integrals C 2 {t) and C${t). 
However, we are interested only in the long-time behaviour of P(t), and can therefore 
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take advantage of tools from asymptotic theory of Laplace and Fourier transforms 
[32J [33], to obtain the asymptotic form of P(t). The functions C 2 (i) and C 3 (t) can 
be written as Fourier integrals of the type Cj(t) = J °° dEfj(E)e~ tEt , with j = 2, 3. By 
considering the explicit form of the eigenf unctions ipo, 4>e, one gets for small values of 
E: 

I 2 ^E^ 2 , I 3 ^E^ 2 E' P/ \ (7) 

with j3 = (a — 3)/2. Hereafter, we will only consider initial conditions ^o(x) for which 
ao 7^ and the spectral decomposition involves continuum energies that extend down 
to E = 0. 

For initial states that take non-vanishing values only in the well [-L, L], ag ~ E 13 
in the limit E — > [27]. Equation ([7]) furthermore allows to infer / 2 (P) ~ K(a) E 13 
and fz{E) H(a)E a ~ 2 , for small values of £7, where and if (a) are prefactors 

that depend also on the initial condition. If the initial state ^o(x) has nonvanishing 
values outside the well region [— L, L], the above results still hold, provided that, for 
large \x\ ^> L, ^o(x) exhibits power-law ~ \x\~ a , Gaussian ~ exp(— ax 2 ), or stretched 
exponential ~ exp(— ax b ) decay. The leading contribution to the asymptotic long-time 
behaviour of the survival probability is then found to be associated with the C 2 (t) term 
above, and reads: 

P(t) pa |a | 2 /i + 2a P 2 t(a _ 1 1)/2 , ast^oo, (8) 
P 2 « TTT^TT 2 cos (I " !)) • (9) 



(2tt)/ 3 + 1 V2 

Through the dependence of K(a) on the spectral expansion of ^o(x), via / 2 or C 2 , the 
asymptotic decay of P(t) also depends on the initial condition. Thus, as anticipated in 
the introduction, we have shown that the potential of equation ([1]) induces an algebraic 
decay in the survival probability P{t), for generic, confined initial states. In contrast, for 
an initial state ^o(x) with = for E < E c (for some E c > 0), one can show that the 
survival probability has an exponential cut-off e~ Ect at large t, which is a manifestation 
of the power-law decay arising from continuous spectral components arbitrarily close to 
the ground state energy. 

Figure [1] shows P 2 for the initial condition ^o(x) = 5(x). For some values of a, e.g. 
in a = 1 and for a = 3, P 2 vanishes. The oscillatory behaviour of P 2 stems from the 
trigonometric functions in ([9]), while the fact that P 2 becomes negligible for a > 7 is 
induced by K(a), and is therefore related to the initial conditions as used in figure [TJ 
For vanishing P 2 , one has to consider the next leading time dependent term in the 
asymptotic expansion of C 2 and C 3 . 

In reference [23J, according to equation (35) therein, a specific initial state was 
chosen whose expansion coefficients a# show a power-law behaviour, although with 
an exponent different from ours, for E — > 0. This implies a power-law for the survival 
probability which is different from our result of equation (JH1). Moreover, for the potential 
well considered in [23] there does not exist a bound state at E = 0, leading to a decay 
without lower bound, or in other words without control, in contrast to our equation (JE|) . 
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Figure 1. Prcfactor P2 of the algebraically decaying term in (7), for the potential (JXJ) 
and different initial conditions: (i) Dirac delta &o(x) — S(x) (solid red line); (ii) water 
bag ^o(x) = l/(2L) in [-L,L] and ^o(x) — outside (dashed blue line); (hi) Gaussian 
distribution ^o{x) = exp(—x 2 /(2s))/iy2Trs, with widths s = 0.1 (dotted green line) 
and (iv) s = 1.0 (dash-dotted magenta line). 



2.2. The role of the potential tails 

In this section we will show that the mere absence of a gap between the null eigenvalue 
and the continuum part of the spectrum, is a necessary but not sufficient condition in 
order to observe a power-law decaying survival probability. 

In general, any Schrodinger potential that asymptotically decays to zero like 
would have a continuum part of the spectrum attached to the null eigenvalue. Let us 
now therefore consider the more general case of a Schrodinger potential 

f -V if |ac| < L, 
V^x) = I /i < 2 (10) 

[ Vx/lx^ if \x\ > L, 

where L, Vq, V\ and fi are real positive constants. Let us call {ipo^, ip^} the 
eigenfunctions associated to such potential. In this case, we can only obtain the 
eigenfunction associated to the null eigenvalue: 

%l,f(x) = Aoy/xKs^d^x 1 -^ \x\ > L 



The eigenfunctions relative to the continuum part of the spectrum are not known. The 
parameters L, V , V\ and \i can be chosen in such a way that the spectrum contains 
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one single discrete eigenvalue E = and a continuous part for E > 0. As a result, 
the parameters L, V , V\ and \i are not independent. In fact, the continuity of d x ip in 
x = L provides a relation between them 

L 1+ ^ 2 Ks^^L 1 -^ (l + 2L v /Vbtan(Lv / H)) + 

- L 3/2 y/ViK s+1 ^V^L 1 -^ (12) 

- L 3/2 v^^-i (s^ViL 1 -^ 2 } = 

In the following we consider L, Vq and fi as independent parameters and will obtain V\ 
by numerically solving the above equation. 

We will prove below that the survival probability of the above process of equation 
( fit)]) is not in general power-law decaying. In fact, we will show that the C<i term may 
eventually decay like a power-law for large time values only for very specific initial 
conditions. Let us make two preliminary observations. First, it is worth mentioning 
that by using the orthogonality relation: 

/>oo 

/ dE^ ) (x)^f(x') = 8(x-x') (13) 
Jo 

one can prove the following identity: 

dEI 2 (E)^ E (x) = ^\x), (14) 



For I2 as defined in equation (J6J). Secondly, we notice that in order to have C^r) « 
l/r a+1 one must have that a E 12(E) « E a for small energy values. 

Let us now consider the large x behaviour of J °° dEl2(E)ipE(x)- For large values of 
x the potential vanishes thus ip^ (x) ~ E^ 1 ^ 4 e lyr ^ x . Therefore, by using equation ffl4l 
and equation (12.01) in chapter 3 of reference [32J, the ansatz 12(E) ~ E a would give: 

dE h(E)Mx) ^ ^j- 2 ( 15 ) 

This result implies that ipQ (x) should decay according to a power-law, given equation 
( |T4l) . However, by using equation ( TTTT) one has that: 



^\x) w x^ /4 exp(-^^x 1 -^ 2 ) (16) 

showing that 12(E) cannot behave as a power-law for small values of E when fi < 2. It is 
worth mentioning that when /i = 2, then the above ansatz holds true with a = (a — 5)/4. 
In this case, one would get /(x) ~ x 1 ""/ 2 , that is in agreement with the fact that 
ipo (x) = x~ a l 2 for large x values. 

We have therefore shown that 12(E) is not growing like a power-law for small energy 
values. Therefore in order to have a power-law decay in C*2(t) one has to engineer 
appropriate initial conditions such that a# 12(E) behaves like a power-law for small 
energy values. In conclusion, the absence of an upper bound for the time-scale is a 
necessary but not sufficient condition in order to observe a power-law decaying auto- 
correlation function. This also implies that the decay of the survival probability is faster 
for fi < 2 than for the case with [i = 2. 
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Figure 2. Sketch of the smooth potential V sm (x) for a — 2.8,7 — 0-68 (black solid 
line), and of the initial state ^o(x) (red dashed line), equation ([!§)) . The thick blue 
line shows the ground state energy level within the well. 



3. Numerical results 

3.1. Numerical Simulations of the model 

Because of its discontinuities, it is difficult to reproduce Vs(x) experimentally or in 
numerical simulations. We therefore consider the smoothed version V sm (x) defined as [26] 

4 [x A + 

in our subsequent numerical tests. Also this potential has a continuous spectrum 
attached to the null eigenvalue, and decays as x~ 2 at \x\ — > 00. We are interested in the 
relaxation properties of wave-packets prepared inside the potential well (see figure [5]). 
We define the survival probability over the region bounded by the potential maxima 
at x c : P(t) = f* c Xc dx\^/(x,t)\ 2 , with x c = 7 ((a + 6)/(a + 2)) 1 / 2 . Our observable is 
the approach of the survival probability to its asymptotic constant value P^ = |a | 2 /i 
described by equation ([8]), i.e., the quantity 

P s {t) = \P{t)-P 00 \. (18) 

Even if we simulate the evolution of the initial wave packet only in one dimension, the 
numerical computations take quite a long time for several reasons. First of all, we need 
to propagate sufficiently far out into the tails of the potential, while - at the same time 
- the expansion of the decaying parts is rather fast due to the high energy components 
of the initial state. This is also the reason why simple absorption methods at the 
numerical boundary do not work very well since those are typically adapted to absorb 
just a small window of energies with sufficient precision. Moreover, in order to estimate 
the asymptotic probability P^, which is not analytically known for the smooth potential 
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Figure 3. Survival probability Ps{t), for a water-bag like initial state of equation (fT9|) 
with a = 100 and 7 = 0.68, and increasing values of a = 2.8 (black dashed line), 3.2 
(blue dot-dashed), 4 (thick green solid line). The straight black line fits are chosen 
according to the prediction with exponent (a — l)/2 from equation (7). 



V sm , we must propagate considerably longer than shown in the following figures, for 
which we approximate by P(t max ), with t max = 140 . . . 200. 

For the actual observation of the predicted asymptotic behaviour the precise form 
of the initial state prepared in the well is not crucial, provided it is given by a coherent 
superposition of energy eigenstates from a continuous energy range including the ground 
state energy We consider a smooth initial state 

»„(*) = , - » - , (19) 

[1 + exp ((x - x c )a)\ [1 + exp (—(x + x c )o)\ 

where N is a normalization constant. The initial state ( I19p . shown in figure I2J is 
numerically propagated in real time using an implicit norm-preserving Crank-Nicolson 
integration scheme [34J, which also controls the boundary conditions very well [35J. 

Figure [3] reports the clear algebraic decay as predicted by equation (7), for various 
values of a which induce different algebraic decay exponents. The power-law decay 
emerges the earlier, the larger the weight factor P 2 in (7). The asymptotic value is 
approached from above (a = 2.8) or from below (a = 3.2,4), depending on this factor's 
sign (compare figure 1). In the latter case, a change of sign of P(t) — at finite times 
induces the discontinuities in the first derivative on the logarithmic scale of figure |3j 
For fixed a, the power- law scaling regime may be further enhanced by adapting the 
precise form of the initial state (e.g. by the above parameter a), and also by changing 
the second parameter 7 of the potential of equation ( ITT]) . Both enter our approximate 
formula for the prefactor P 2 , see equation (jHJ) and figured) 
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3.2. The role of the tail in the potential 

In the following we will show that the predicted power-law decay is robust with respect 
to realistic experimental situations. The only parameter which needs to be controlled 
very well is the exponent of the potential tails, i.e. /i in the asymptotics of the potential 
scaling as ~ x~^ . We use the form 

V ( x )-^ x2(1 + a/2) ~ 72 (20) 

»W W - 4 ( x 2 + 7 2)^/2+l ' 1 ' 

which for \i = 2 recovers V sm (x) from equation ( TT71) . 

For /i ^ 2, the dynamics of the system will be qualitatively very different, as we 
show in figure H] directly for the survival probability. Instead of an asymptotically slow 
saturation toward the ground state within the potential, we observe a fast decay - with 
approximately constant slopes - for all values fi < 2 along the same time scales as in 
figure EJ 

As described above around equation (fl8"l) . it is computationally hard to estimate 
the saturation value of the survival probabilities, since the ground states are not known 
analytically and the data for /i ^ 2 implies that a reliable numerical estimate will be 
possible only after much longer propagations than shown in the figures. This may cause 
the wrong impression suggested by the main panels of figure HJ the ground state is 
approached in shorter absolute time for \i = 2 as compared to the other values of fj,, yet 
what counts is the rate of approach which is always smaller for \i — 2. This is what we 
intend as 'slower' decay in the latter case and the reason for plotting the slopes defined 
by 

(P(120)-P(t))/At 

(P(120) - P(20))/100 ' 1 ' 

in the insets, with At = 120 — t. In order to obtained comparable values, we divided 
by the extremal values, see the denominator of (I2ip . which all occur at the minimal 
time t = 20 (chosen after non- universal transients at still smaller times). Only in the 
special case of /i = 2, the relaxation to the steady-state value becomes much slower with 
time, see the fast decreasing slopes in the insets (black solid line). The shown numerical 
results confirm the theoretical discussion above in section I2.2[ predicting in essence a 
faster approach to the ground or equilibrium state for /i < 2. Consequently, a possible 
experiment should control the exponent to be equal to two at least at 2-3 significant 
digits, in order to observe the power-law decay discussed in sections 12.11 and 13.11 



3.3. Experimental realization of the model 

Let us conclude with an experimental protocol to test our prediction (jSJ). We have in 
mind a Bose-Einstein condensate prepared in an optical trap which then is exposed 
to a potential of the form of equation ( !T7l) while the trap is relaxed. The potential 
may hereby be created optically, for instance, with a fast moving laser beam [18] or by 
holographic techniques [36J. 




Figure 4. Survival probability P(t), for a water-bag like initial state of equation (fT9|) 
with a = 60, for well parameters 7 = 0.68, a = 2.8 (a) and a = 4 (b), and asymptotics 
with exponents fi = 1.5 (green dotted line), /i = 1.9 (blue dashed line) and /i = 2 
(black solid line), respectively. The insets show, for t = 20 . . . 120, the decay rates or 
slopes defined in @IJ): (P(120) - P(t))/At vs. 100 - Af = t - 20, with A< = 120 - t. 
Hence, we observe a relaxation with a fast decreasing slop as At — > 0, in the special 
case of /i = 2. In all the other cases the slope remains almost constant over the shown 
times, implying a much faster decay with time. 



A Gaussian initial state of the form ^o(x) = N exp(— x 2 /(4of rap j)) (N a 
normalization constant) is initially prepared in a harmonic trap with characteristic 
oscillator length Otrap.i- Then the smooth potential (ITT)) is switched on, while we 
switch off or relax the trap abruptly to a shallow confinement characterized by the 
harmonic oscillator length cr tra p,f 3> cr traPii . In the case of figure [5] we consider values 
from cr tra p,f = 100 to er trap ,f = 800, and o" traPi i = 0.5. In the absence of a trap (thick 
blue lines), the Gaussian initial state exhibits a behaviour similar to that of the smooth 
water-bag state from above, equation (fl9|) . A shallow trap with large cr t ra P ,f manifests at 
long times, by an exponential cut-off of Ps{t), as shown in figure El since the potential's 
asymptotics are changed by the trap. This induces a spectral gap between the ground 
state energy and the continuum component, while it is specifically the absence of the 
gap which is responsible for the algebraic decay, as discussed above. The steeper the 
confining trap potential, the larger the trap-induced spectral gap, and the shorter the 
time interval over which an algebraic decay can be observed, before the exponential 
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Figure 5. (a) Survival probability for a confining potential with a — 4 (i.e. 
x c = 0.88), and otherwise the same parameters as in figure |31 The particle is here 
prepared in a Gaussian initial state with a spatial width <7t ra p,i = 0.5, and evolves in 
the potential (|17p without further perturbation (thick blue lines), or in the presence 
of an additional harmonic trap of width <7t rap ,f = 100 (crosses), 150 (red dashed), 200 
(green dot-dashed), 400 (dotted), and 800 (black squares). Inset: long-time behaviour, 
(b) Data from (a) shown in the form of Ps(t), equation (fT8|) . for direct comparison 
with figure [3] and the power law of equation (7), with (a — l)/2 = 1.5 (thin solid line). 



cut-off. Indeed, we observe such a continuous degradation of the asymptotic law with 
exponent (a — l)/2 = 1.5 when making the additional confinement steeper. However, 
the power-law trend is still clearly visible, over at least one order of magnitude, even in 
the presence of the steepest trap (with cr traPi f = 100 in figure |5]). 

One may ask about the role of interactions between the weakly interacting atoms 
of a Bose condensate. Those effectively scale in a mean-field approximation with the 
number of atoms, which may be controlled and possibly be reduced [H]. We modelled 
the evolution of an initially well confined condensate, following the above protocol, using 
a one-dimensional Gross-Pitaevskii equation. In this approach the interactions are taken 
into account by a nonlinear density- dependent term in the Hamiltonian [HIUS, 35J. A 
weak repulsive nonlinearity may actually stabilize the evolution and - to some extent 
- reverse the effect of a weak confinement during the relaxation (c.f. [37] for a similar 
effect). This can be seen in our final figure [61 where the prefactor of the nonlinear term 
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Figure 6. Direct and shifted survival probabilities Pit) (a) and Ps(t) (b), 
respectively, taken from figure [5] without any trap (thick black solid line) and for 
otrap.f = 100 (crosses), 150 (thin red dashed), 200 (thin green dot-dashed) for 5 = 0, 
in comparison with the same cases and a weak repulsive nonlinearity g = 5xl0~ 3 and 
no trap (blue dotted) or <7t ra p,f = 100 (crosses connected by thin line), 150 (thick blue 
dashed), 200 (thick blue dot-dashed). 



is denoted by g in our dimensionless unitq§|. 
4. Conclusions 

In summary, we propose a setup for which one may observe a power-law decay of the 
asymptotic survival probability Ps(t) in a controlled manner. Such an anomalous decay 
is readily realized by preparing a quantum wavepacket with a sufficiently broad energy 
range around the ground state energy E = in the potential wells given by Eqs. ([I]) and 
( fl7l) . The slow algebraic relaxation toward this equilibrium state arises from the spectral 
properties and the population of eigenmodes of the system. The latter are encoded in 
the coefficient C 2 defined in equation (JSJ). Hence, the predicted behaviour is of purely 

§ Our dimensionless potentials are given in units of energy /m 2 . The three-dimensional atom- 
atom interaction strength can be reduced to an effective one-dimensional parameter, provided a 
strong transverse/radial confinement is experimentally achieved |14j . This latter parameter guy = 
2hu> Ia ,dQ'sN a , with the scattering length as, the radial confinement frequency w r ad and the number of 
atoms in the Bose condensate N a , can be expressed without dimensions using, e.g., just the single scale 
given by the radial confinement: Eq — hu! la d and xq = -\/2/?kj ra dM, where M is the single atom mass. 
This gives g — gm/iEoXo) = 2N a as/xo- Here the wave function is normalized to one as in our single 
particle computations otherwise used in the paper. 
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quantum origin (not induced, however, by dynamical or Anderson localization jUE] 
as discussed in j38j, nor by semiclassical arguments [39J). The exceptional control of 
state-of-the-art experiments with ultracold atoms [T5llT8l42T| l36] offers the possibility to 
observe our predictions, following, e.g., our protocol for obtaining the data of figures [5] 
and El 
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